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ABSTRACT 

Context. In the recent literature there has been some doubt as to the reliability of CO multi-transitional line observations as a mass- 
loss-rate estimator for AGB stars. 

Aims. Using new well calibrated CO radio line observations, the main aim of the work presented here is to carefully evaluate the 
reliability of CO mass-loss-rate estimates for intermediate- to high-mass-loss-rate AGB stars with different photospheric chemistries. 
Methods. Mass-loss rates for 10 intermediate- to high-mass-loss-rate AGB stars are derived using a detailed non-LTE, non-local 
radiative transfer code based on the Monte-Carlo method to model the CO radio line intensities. The circumstellar envelopes are 
assumed to be spherically symmetric and formed by constant mass-loss rates. The energy balance is solved self-consistently and the 
effects of dust on the radiation field and thermal balance are included. An independent estimate of the mass-loss rate is also obtained 
from the combination of dust radiative transfer modelling with a dynamical model of the gas and dust particles. 
Results. We find that the CO radio line intensities and shapes are successfully reproduced for the majority of our objects assuming 
a constant mass-loss rate. Moreover, the CO line intensities are only weakly dependent on the adopted micro-turbulent velocity, in 
contrast to recent claims in the literature. The two methods used in the present work to derive mass-loss-rates are consistent within a 
factor of ~ 3 for intermediate- to high-mass-loss-rate objects, indicating that this is a lower limit to the uncertainty in present mass- 
loss-rate estimates. We find a tentative trend with chemistry. Mass-loss rates from the dust/dynamical model are systematically higher 
than those from the CO model for the carbon stars and vice versa for the M-type stars. This could be ascribed to a discrepancy in the 
adopted CO/H^-abundance ratio, but we caution that the sample is small and systematic errors cannot be excluded. 
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of AGB stars are, at least partly, dust-driven, i.e., the mass-loss 
from the star is driven by radiation pressure on grains and the gas 
is dragged away from the s tar by momentum-coupling between 
the grains and the gas (e.g jLamers & Cassinellil[l999l) . The ef- 
fects of the gas-dust interaction on the mass loss and wind for- 
mation have been studied by several authors (e. g.. iHabindl 1 9961 : 
ISimis et aDl200lh . Time-depend ent dynamical models studying 
the effects of stellar pulsations dBowen|[l988l) . a nd further in- 
cludin g detailed time-dependent dust formation dHofner et al.l 
119951) . have progressed to dynamical models in which the mo- 
tion of the gas and du st particles is calculated and the radia- 
tive transfer is solved (ISandin&Hofnerl 120031) . These models 
work well for carbon stars, and support the dust-driven wind sce- 
nario, but the formation of silicate grains close enough to the star 
to drive mass-loss rates of t he observed magnitude still poses 
a pro blem for M-type stars dWoitkd [20061 : iHofner & Andersen] 
12007b . 

Most likely, the averag e mass-loss rate increases as the 
star e volves along the AGB (Habing 1996; Sch oier & Olofssonl 
1200 lh . In addition to this steady increase, mod ulations of th e 
mass-loss rate on shorter time scales are present (Habing 1996). 
Indications of this can be seen in observations of circumstel- 
lar CO line emission arou nd some AGB star s with double- 
comp onent line profiles dKnapp et alj 119981: IWinters et al.l 
120031) . and in observations of stellar and ambient Galactic 
light scattered by dust i n circumste llar envelopes (CSEs) 
dMauron & Hugginsll2000l 1200 ll 120061) . A process which can 



1. Introduction 

Mass loss is the single most important process during the 
final evolution of low- and inte rmediate-mass s tars on the 
asymptotic giant branch (AGB) dBloeckerl I1995I) . The mass 
loss determines the lifetime on the AGB, a nd thus the max- 
imum luminosity and radius reached (e.g., iMarigo & Girardil 
|2007|) . Furthermore, since AGB stars are major contribu- 
tors of mass to the inter stellar medium dSedlmavrl Il994t 
iDorschner & Henning||l995l) . and to the integral luminosity of 
intermediate-age stellar systems, an understanding of the mass- 
loss phenom enon is also important for studying extra-galactic 
populations dBattinelli & Demersl 2004I) and for the chemical 



evolution of galaxie s I Schroder & Sedlmavrl d200ll) : see also 
iKarakas & Lattanziol d2007l) for a recent paper on yields from 
AGB stars] . T he existence of mass loss is well established 
dWillsonll2000l) . but its details are less well known. This is un- 
fortunate since even a modest mass-loss-rate change, e.g. by a 
factor of two, will have a profound effect on the evolution of the 
star, its nucleosynthesis, and its mass return to the interstellar 
medium dForestini & Charbonn el 1997). 

Although mass loss on the AGB can reach rates as high as 
lO~ 4 M0yr _1 , the wind speed is quite low (<30kms _1 ). The 
temperature of the outer parts of the atmosphere is low enough 
for dust to condense. These facts, combined with a high stellar 
luminosity (~ 10 4 L ), have lead to the conclusion that the winds 
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have a strong effect on the mass-loss properties of an AGB-star is 
the He-shell flash (or thermal pulse) (lOlofsson et al II 19931 l2000t 
ISteffen & Schonbernerll2000t ISchoieret al.ll2005al) . Other pro- 
cesses may be important during the very last stages when the star 
sheds the remaining few tenths of a solar mass of its hydrogen 
envelope (during a f ew hundred years) before it leaves the AGB 
dSkinner et al.ll 1997b . Finally, the presence of a binary compan- 
ion may influence the magnitude as well as the characteristics 
(e.g., the geometry) of the mass loss. 

To make progress in the study of AGB stars, it is important 
to establish reliable mass-loss-rate estimators, which preferably 
are useable over a wide range of mass-loss rates and time scales, 
and for different chemical types (M-, S-, and C-stars). Estimating 
AGB mass-loss rates through observations of circumstellar CO 
radio line emission, in combination with modelling of the radia- 
tive transfer, is considered to be one of the most reliable meth- 
ods. CO line emission has also been used in a number of stud- 



Table 1. The sample with spectral types, variable types, periods 
(P), luminosities (L+), and distances (D). 



ies aimed at tracin 


g the mass-loss history (Schoieret alj|2002l 


iKemper et alj|2003 


jTevssier et al.ll2006tlDecin et al.ll2007l). and 



at ide ntifying drastic mass-loss-rate changes at the end of the 
AGB dHeskeetal .119 90; Justt anont et al.ll996l). Detailed studies 



of the mass-loss propertie s of carbon stars dSchoier & Olofsson 



2001; ISchoieret all l2002l), M-type A GB stars (lOlofsson et al 
2002t IGonzalez Delgado et al l [2003b . and S-type AGB stars 



dRamstedt et al.l 120061) have been performed on samples domi- 
nated by low-mass-loss-rate objects. It was concluded that for 
this type of objects the CO radio line emission is a good mass- 
loss-rate estimator, although not entirely free of problems. 

It has been known for some time that for the higher mass- 
loss-rate objects the situation is less clear {Sahai 19961 iKastnerl 
fl99l ISchoier & Olofssonll200lh . A higher mass-loss rate leads 
to a decrease in the drift velocity between the dust and the gas 
and therefore a less efficient heating of the gas through colli- 
sions with dust grains. Since CO line emission is also a dom- 
inant coolant, the effect of an increased mass-loss rate is in- 
creased cooling of the circumstellar gas. This, in combination 
with saturation of the spectral lines, leads to a situation where 
the integrated CO line intensities are less dependent on the mass- 
loss rate when it exceeds about 10~ 5 M yr~' ; the higher the fre- 
quency of the lines the lower the mass-loss rate when the ef- 
fect sets in. This makes CO line emission a questionable mass- 
loss-rate estimator at mass-loss rates in excess of 10~ 5 M yr _1 , 
an unfortunate effect since we know that mass-loss rates up to 
10^ 4 M Q yr _1 are reached during the very final stage, and this 
is also when the stars make their largest contribution to the re- 
turn of n uclear-processed mate rial to the interstellar medium. In 
addition, IKemper et al.l d2003l) concluded that standard circum- 
stellar radiative transfer models could not explain some CO line 
intensity ratios, and that the mass-loss-rate estimate is sensitively 
dependent on the (unknown) local turbulent line width. 

An alternative way to estimate the mass-loss rate, not subject 
to the same problems, is to combine modelling of the dust radia- 
tive transfer with solving th e coupled equations of motion for the 
gas and the dust particles dHabing et al.lll994b . The gas expan- 
sion velocity is a strong function of the dust-to-gas mass-loss- 
rate ratio and therefore a solution of the kinematical equations 
that fits the gas expansion velocity and the dust optical depth, 
obtained through dust radiative transfer modelling, results in an 
independent estimate of the mass-loss rate. 

Presented here are a detailed radiative transfer modelling 
of primarily new, well-calibrated, and high signal-to-noise ratio 
observations of CO / = 1 -> 0, 2-* I, 3 — > 2, 4 — > 3, and 6 — > 5 
emission from a sample of 10 well-studied high-mass-loss-rate 
AGB stars, and the resulting mass-loss-rate estimates. In addi- 
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Mira 


500 


10400 


300 


IRC-10529 


M 


Mira 


680 


15100 


890 



tion, dust emission modelling combined with a dynamical model 
give independent estimates of the mass-loss rates. The dust emis- 
sion modelling also provides the dust information necessary for 
the CO modelling, e.g., the dust radiation field. By obtaining 
the mass-loss rates for the stars using two independent methods, 
we can evaluate the reliability of the mass-loss-rate estimators, 
and also look for trends with e.g. mass-loss rate and chemistry. 
A similar investigation comparing mass-loss-rates derived from 
ISO-SWS observations of the silicate dust features with mass- 
loss rates from CO modelling found in the literature on a sample 
of low-mass-loss-rate M-type AGB stars has been performed by 
lHeras&Hoiivld2005l) . 

The sample and the observational data are presented in 
Sects [2] and [3] In Sect. |4]the adopted circumstellar model is dis- 
cussed. In Sect. [5] the CO radiative transfer model is described 
and CO as a mass-loss rate estimator is discussed. In Sect.|6l we 
describe how dust radiative transfer and the dynamical model 
is used to estimate mass-loss rates. The results are presented in 
Sect. [7] and outlined for each individual object in Sect. [8] In 
Sect.|9]we discuss and compare the two methods, and in Sect.lTOl 
we give our conclusions. 



2. The sample 

The sample, presented in Table [1] consists of 10 well-studied 
AGB stars, with intermediate to high mass-loss rates. They are 
chosen with a spread in mass-loss rate, to allow studies of possi- 
ble trends. Dependences on chemistry can also be investigated, 
since half of the stars were chosen to have a photospheric C/O- 
ratio >1 (carbon stars) and the other half C/0<1 (M-type). 

Bolometric lumin osities are e stimat ed using the period- 
luminosity rekdtion_ofjFeast et al. (2006) for the carbon stars, 
and that of lWhitelock et alj d 19941) for the M-type stars. Periods 
are mainly tak en from the Comb ined General Catalogue of 
Variable stars dSamus et al.l 120041). The periods for LP And , 
AFGL 3068, and IK Tau are from iMen'shchikov et alj (I2006T) . 
iLe Bertrd d 19921) . and IGonzalez Delgado et all d2003l) . respec- 
tively. Distances are derived by fitting the spectral energy dis- 
tribution (SED) calculated with DUSTY (Sect. |6J} to ob- 
served fluxes from the near-IR to the mm-range (Sect. 13.21 and 
AppendixlBt. 
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Table 2. Main-beam efficiencies (77 m b), and main-beam 
FWHM:s (0 m b), at the frequencies of the CO lines. 
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0.43 


33 


oso 


J = 2 -> 1 


230.538 


0.66 


21 


JCMT 


J = 3 -> 2 


345.796 


0.62 


14 


JCMT 


7 = 4^3 


461.041 


0.5 


11 


JCMT 


J = 6 -> 5 


691.473 


0.3 


8 


JCMT 



Table 3. Integrated CO(7 — > J — 1) intensities, J T m \,dv, in 
Kkms 1 . 
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35 


63 


40 


79 




AFGL 3068 


48 


72 


73 


93 


55 
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20 


61 


71 


149 
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31 


61 


80 


79 
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44 


81 


70 


84 


23 
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46 


95 


125" 


130 




IRC-10529 


16 


45 


27 


22 





" JCMT archival data 



3. Observational data 

3.1. New observations of CO radio line emission 

When estimating the physical properties of the CSEs, e.g., mass- 
loss rates and kinetic temperatures, it is very important to hav e 
well-calibrated data for all lines dSchoier & Olofssonl |2001|) . 
CO /= 1 — >0 (115.271 GHz) data was obtained in December 
2003 using the Onsala Space Observatory (OSO) 20 m tele- 
scoped The James Clerk Maxwell TelescopeQ (JCMT) located 
on Mauna Kea was used to collect CO 7 = 2— >1 (230.538 
GHz), 7 = 3 — > 2 (345.796 GHz), / = 4-* 3 (461.041 GHz), and 
7 = 6 — > 5 (691.473 GHz) line emission from the sample sources. 
The observations were performed in June, July, and October 
2003. 

All observations were made in a dual-beamswitch mode, 
where the source is alternately placed in the signal and the refer- 
ence beam, using a beam throw of about 2' at JCMT and 1 1' at 
OSO. This method produces very flat baselines. Regular point- 
ing checks made on strong continuum sources (JCMT) or SiO 
masers (OSO) were found to be consistent with the pointing 
model within «3". 



1 The Onsala 20 m telescope is operated by the Swedish National 
Facility for Radio Astronomy, Onsala Space observatory at Chalmers 
University of technology 

2 The JCMT is operated by the Joint Astronomy Centre in Hilo, 
Hawaii on behalf of the present organizations: the Particle Physics and 
Astronomy Research Council in the United Kingdom, the National 
Research Council of Canada and the Netherlands Organization for 
Scientific Research. 



The data was reduced in a standard way using XS0by remov- 
ing a first order polynomial fitted to the emission-free channel 
ranges, and then binned in order to improve the signal-to-noise 
ratio. The intensity scale is given in main-beam-brightness tem- 
perature, r m b = r£ /i] m b, where T* is the antenna temperature 
corrected for atmospheric attenuation using the chopper-wheel 
method, and Tj m b the main-beam efficiency. The adopted main- 
beam efficiencies are given in Table [2] The uncertainty in the 
absolute intensity scale is estimated to be about +20% for all 
lines except the CO 7 = 6 — > 5 lines, where we estimated a higher 
uncertainty of about +30%. 

The observed spectra are presented in Figs [3] |4] and IC.ll - 
IC.8I and velocity-integrated intensities are reported in Table [3] 

3.2. Dust continuum emission 

The SEDs are constructed from JHKLM-band, sub-millimetre, 
and millimetre photometric flux densities from the literature, 
together with IRAS data. SCUBA data (at 450 and 850 /an) 
was downloaded from the JCMT archive for CW Leo, AFGL 
3068, and RW LMi. All fluxes and references are presented in 
AppendixlBl 

To check the dust radiative transfer model results, IRAS 
and ISO spectroscopic data were used. IRAS-LRS spectra were 
downloaded^ for all stars in the sam ple. An absolute calibr ation 
correction was app lied according to IVolk & Cohen! d!989l) and 
ICohen et alj (199% . The ISO LWS and SWS Highly Processed 
Data Products (HPDP) were downloaded and not processed fur- 
ther. LWS spectra were available for all stars, except GX Mon. 
SWS spectra were found for TX Cam and WX Psc (SWS spectra 
for carbon stars were not searched for). For IK Tau, only a grat- 
ing scan covering the wavelength range 29^-9 fim was avail- 
able, and for IRC-10529, the SWS spectrum is very noisy above 
28 /urn. 



4. The circumstellar model 

The CSE is assumed to be spherically symmetric and formed 
by a constant mass-loss rate. It is assumed to be expanding at a 
constant velocity, derived from fitting the CO line widths, and to 
have a micro-turbulent velocity distribution with a doppler width 
of 0.5 km s (e-folding value) throughout. In addition, a thermal 
contribution to the local line width is added, based on the derived 
kinetic temperature of the gas. The density structure is obtained 
from the conservation of mass. 

Even in this simple model there are a number of uncertain- 
ties that will affect the strengths and the shapes of the CO lines 
and the shape of the SED, e.g., the thermodynamics of the gas, 
the CO chemistry, the dust composition, the optical properties 
of the dust, and the grain size distribution. In addition, there can 
be major deviations from this model due to e.g. non-isotropic, 
time-variable mass loss, and a clumpy medium. Finally, the ob- 
jects are long-period variables and the data are sometimes taken 
at different epochs. In many cases it may be impossible to disen- 
tangle which of these assumptions and effects, if any, are affect- 
ing the emerging radiation. This taken into account, the most rea- 
sonable approach is to choose the simplified circumstellar model 
and analyze all stars using the same method. 



3 XS is a package developed by P. Bergman to reduce and analyze 
a large number of single-dish spectra. It is publicly available from 
ftp : //yggdrasil . oso . Chalmers . se 

4 http://www.iras.ucalgary.ca/~volk/getlrs_plot.html 
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5. CO line modelling 

Only a short description of the radiative transfer analysis of the 
circumstellar CO radio line emission is given here, since the 
adopted model has been de scribed in detail in p revious articles 
Schoier & Olo fsson (2001) for carbon stars, and lOlofsson et al.l 



20021) for M- type starsl . 



The non-LTE radiative transfer code is based on the Monte 
Carlo method. 41 rotational levels are included in each of the two 
lowest vibrational states (v=0 and v=l) in the excitation anal- 
ysis of the CO molecule. Energy levels and radiative transition 
probabilities, as well as collisional r ate coefficients for coll isions 
between CO and H2, are taken from lSchoier et al.l (l200 5bFl An 
ortho-to-para ratio of 3 is adopted when weighting together colli- 
sional rate coefficients for CO in collisions with ortho- and para- 
H2. Radiation from the central source (assumed to be a black- 
body), thermal dust grains (Sect. lOT i. and the cosmic microwave 
background is included. 

The energy balance equation for the gas is solved self- 
consistently, by including line cooling from CO and H2 and cool- 
ing due to the adiabatic expansion of the gas. The dominant gas- 
heating mechanism is collisions between the H2 molecules and 
the dust grains. Photoelectric heating is also included and it is 
mostly important in the outer parts of the CSE. The free parame- 
ters describing the dust are combined in the parameter, h, defined 
as 



h = 



¥ \/2.0gcirr 3 \/0.05/mi 



0.01 



Pg 



(1) 



where *P is the dust-to-gas mass-loss-rate ratio and p g and 
fl g are the average density a nd radius of an ind i vidual dust 
grain, respectively. Following Sch6i er & Olofssonl d2001l) . we 
have adopted an average efficiency factor for momentum trans- 
fer (see Sect. 16.21 ), which is constant throughout the CSE, 
(Grp) 1 =3 X 10 -2 . However, the particular choice of (<2rp) is not 
so critical. As long as the product /i(2rp) ,3 ^ 2 is kept constant (and 
h is a free parameter), the model results are essentially the same 
for intermediate and hi gh mass-loss-rate objects, f or a given M 
(see Eqns (4) and (6) in lSchoier & Olofssonl d200 lb ). 

The radial distribution of CO is estimated using the model 
presented in Mamon et al. (1988). It includes photodissociation, 
shielding, and chemical exchange reactions. The photospheric 
abundance of C O relative to H2, is assume d to be 1 x 10 -3 for 
the carbon sta rs dZuckerman & Dvck|[T9 86) and 2 X 10 -4 for the 
M-type stars dKahane & Jurall 19*9*4 ). The luminosity and the in- 
ner radiu s of the CSE were taken fro m the SED fitting. It was 
shown in Sch oier & Olofssonl d2001l) that, in the intermediate- 
to high-mass-loss-rate stars, the CO radio line intensities are es- 
sentially unaffected by a change in these parameters. The main 
reason is that in the higher-mass-loss-rate CSEs, CO is collision- 
ally excited out to large radii. The mass-loss rate and h are the 
remaining free parameters in the CO line modelling. 

6. Dust radiative transfer and dynamical model 

6. 1 . Dust emission modelling 

Dust radiative transfer through a spherically symmetric enve- 
lope possesses scaling properties dlvezic & Elitzur|[l997l) . This 
fact is put to use in the publicly available radiative transfer code 
DUSTY0 adopted here to model the observed continuum emis- 
sion. The scaling properties of the solution makes DUSTY ideal 

5 http://www.strw.leidenuniv.nl/~moldata 

6 http://www.pa.uky.edu/~moshe/dusty/ 



for studying a large sample of stars, since only one model grid, 
although large, needs to be calculated for each dust type. 

The most important parameter in the dust radiative transfer 
modelling is the dust optical depth 



p d (r)dr, 



(2) 



where ka is the dust opacity per unit mass, pa the dust mass den- 
sity, and the integration is made from the inner fa) to the outer 
fa) radius of the CSE. Prompt dust formation is assumed at the 
inner radius. Amorph ous carbon dust grains, with the optical 
constants presented in ISuhl d2000l) . are adopted for the carbon 
stars, and amorphous silicate gra ins are adopted for the M-type 
stars dJusttanont & Tielenslll992l) . For simplicity, the dust grains 
are assumed to be of the same size with a radius of 0. 1 /urn, and 
have a density of 2 g ctrT 3 (carbon grains) or 3 g ctrT 3 (silicate 
grains). 

In the modelling, the dust optical depth, specified at 10 /mi, 
is varied in the range 0.1-10 in steps of -10%. In addition, 
the dust temperature at the inner radius of the dust envelope 
TdfaX an d the stellar effective temperature T e ff, are varied. These 
are the three adjustable parameters in the analysis. The SED 
model result depends only weakly on the other input parame- 
ters, which are fixed at reasonable values. The size of the enve- 
lope is fixed at r e fa=20000 to make sure that the whole enve- 
lo pe extent determined from the CO emission fa; see Table 3 
in lSchoier & Olofssonl d200ll) l is covered. Using an even larger 
size does not affect the result. 

The temperature at the inner radius should be close to (or 
lower than) the dust condensation temperature. In the mod- 
elling it is varied in the range 300-1 500 K, in steps of 100 K. 
The canonical dust condensation temperature of carbonaceous 
dust grains is «1500K, while the condensation temperature for 
the silicate grains is expected to be somewhat lower. However, 
prompt dust formation is an oversimplification. Given the com- 
plex chemistry close to the star, where shocks are present and 
grain growth is important, the dust-to-gas ratio probably de- 
pends on the radius. The dust temperature at the inner radius 
should therefore be treated as a measure of the characteristic 
radius where dust formation is essentially completed. The stel- 
lar effective temperature is varied in the range 1 800-2400 K 
(also in steps of 100 K), characteristic of high-mass-loss-rate 
AGB stars, but lower than that of AGB stars with optically thin 
dust envelopes (e.g., Lambert et al. 1986; Bergeat et al. 2001; 
Heras&Hony2005). 

When scaling the result from DUSTY, the interstellar extinc- 
tion is taken into account. It is calculated fr om the galactic lon- 
gitude for each star, using equation (8) in iGroenewegen et all 
(U992I 



6.2. Dynamical model 

When light is scattered by a dust particle, some of its momentum 
is given to the scattering particle. For particles with an arbitrary 
shape, the efficiency factor for this momentum transfer, depends 
on the orientation of the particle and on the polarization of the 
light. However, for spherical particles, as assumed here, it is in- 
dependent of both. The dimensionless efficiency factor is given 
by the momentum fraction taken from the incident beam, divided 
by the geometrical cross section of a particle: 



Qa = ^^flgPg 



(3) 
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In the CSEs of AGB stars, momentum transfer from photons 
to dust grains pushes the dust outward, and friction between the 
dust and gas particles determines the gas outflow velocity. By 
assuming that the wind is entirely dust-driven, conservation of 
mass and momentum can be used to set up equations for the 
velocity of the dust (v<j) and the gas (v g ) particles, 



dr 



dr 



AncrriA 



GM t 



GM» 



<\itm&v 



(4) 



(5) 



where cr d and m& are the collisional cross section and mass of a 
dust grain, (Q IV ) the radiation pressure efficiency averaged over 
the photon spectrum, L* and M, the stellar luminosity and mass, 
Vdi (=va - v g ) the drift velocity between the dust and the gas, and 
v s the sound speed assumed to be 2 km s _1 dLamers & Cassineilil 
119991). A stellar mass of 1 .5 M is used f or all stars in the sample 
dKahane et aHl2000l:IOlivier et alJl200lh . 

The radiation pressure efficiency Q \ is determined and aver- 
aged over the photon spectrum, 



<&p>to = 



/ Q A B A e^dA 
f B A e~ r *dA 



(6) 



where B A , calculated from the effective temperature of the star, 
is weighted by the absorption determined by the optical depth t a 
derived using the dust opacity and the dust density structure from 
the dust emission modelling (Eqn.|2]i. This includes absorption, 
but not re-emission at longer wavelengths. This simplification 
will probably not have a big effect on the result due to the wave- 
length dependence on the radiation pressure efficiency. DUSTY 
includes an option where the wind structure can be calculated by 
solving the hydrodynamical equations. The drawback compared 
to the dynamical model presented here is that a dust-to-gas ratio 
needs to be assumed in order to get an estimate of the mass-loss 
rate, but as opposed to our model, it includes the re-emission by 
the dust. A comparison between the results from our dynami- 
cal model and the mass-loss-rate estimates derived by DUSTY 
on five of our sources showed that both methods give the same 
result (within ~ 30%) for the same dust-to-gas ratio. 

(Grp) rapidly decreases and reaches its terminal value, 
(2,p)' , at * 20 - 30 R+. Equations <[4]> to (|6]> are solved in an 
iterative way, for different mass-loss rates, until the dust density 
structure is consistent with the optical depth derived in the dust 
radiative transfer model, and a fit to the observed gas expansion 
velocity is found. In this way, an independent estimate of the 
ma ss-loss rate is obtai ned. 

iBains et ail d2003l) investigated maser emission towards four 
AGB stars and derived velocity profiles for their circumstellar 
envelopes. In Fig.Q] our result for IK Tau is overlayed and shown 
to be consistent with their data. 



7. Results 

7.1. CO line modelling 

7.1.1. Strategy for finding the best-fit model 

Only two parameters are adjustable when fitting the model to 
the observed integrated line intensities: the mass-loss rate, M, 
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Fig. 1. The velocity profile for IK Tau from the dynamical model 
overlayed with the velocity data obt ained from mapping of 
maser emission (see IBains et al.l 12003, and references therein). 
H2O masers (circles) and OH maser (squares) have been used. 
The expansion velocity derived from the CO data is marked by 
a diamond. 



and the /i-parameter. Models are calculated for a large number 
of mass-loss rates and /1-parameters for each star and the best-fit 
model is found using^ 2 -statistics. The reduced^ 2 for the best-fit 
model is estimated from 



2 

-^red — 



N ■ 



(7) 



where N is the number of observational constraints, and p the 
number of adjustable parameters (2 in our case). x^ in is the min- 
imum^ 2 



X 



-s 



modi 



^obs.i) 2 



07 



(8) 



where <x, for the high-quality data presented here, is dominated 
by the calibration uncertainties assumed to be 20% for all lines 
except the J=6—>5 lines where they are assumed to be 30% (see 
Sect. ED. 

7.1.2. Mass-loss rates from CO 

We are able to successfully model the CO line emission from all 
stars in our sample, both carbon stars and M-type stars, using 
a constant mass-loss rate. (The J-6— >5 line of WX Psc poses a 
problem that will be further discussed in Sect. l8.2.3l ) The results 
are presented in Table|4] Fig.|2]shows^ 2 -maps for all our sample 
stars. The difficulty in setting an upper limit on the mass-loss 
rate for high-mass-loss-rate stars using CO radio line intensities 
only is apparent. The saturation of the optically thick CO lines 
weakens significantly the dependence of the line intensity on the 
amount of CO in the envelope. 

In some cases, better constraints on the mass-loss rate can be 
set by including the shape of the line profile. Changing the mass- 
loss rate resultsjii_achangjrmphotodissociation radius [see, e.g., 
model by IStanek et al.l (1 1 995l) l . and this may result in a change 
of the line profile. Visual inspection of the line profiles (within 
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Table 4. Mass-loss rates, /i-parameters, and terminal gas expan- 
sion velocities v e , derived by fitting the CO line radiative transfer 
model to the observations. 



Source 


M 


h 


v e 


2 
✓Wed 


N 




[Kr 6 M yr-'] 




[kms- 1 ] 




LP And 


7 


2.0 


14.0 


0.5 


5 


CWLeo 


20 


1.0 


14.5 


0.1 


4 


RWLMi 


5 


1.5 


17.0 


0.7 


4 


V384 Per 


3 


1.0 


14.5 


2.7 


4 


AFGL 3068 


20 


2.0 


14.0 


0.9 


5 


TXCam 


7 


1.0 


18.5 


3.0 


4 


GXMon 


20 


0.5 


18.7 


0.1 


4 


WXPsc 


40 


0.4 


19.3 


1.0 


4 


IK Tau 


10 


0.3 


19.0 


0.4 


4 


IRC-10529 


16 


0.2 


14.0 


3.8 


4 



the range of the innermost contour of the ^ 2 -maps) was used to 
pick the best model presented in Table [4] and when determin- 
ing the error bars in Fig. [7] In principle, it is possible to account 
for the shape of the line profile in a more quantitative way, i.e. 
by dividing the line into several velocity bins and calculating an 
error using ^-statistics. However, this is problematic since the 
resulting ^-values become very sensitive to very small varia- 
tions in the line-profile shape which may not be possible to ac- 
count for by changing the free parameters in the adopted model. 
Also, in cases like TX Cam (Fig.|4|i, a fit to the line-profile shape 
can never be found within the adopted model, either the adopted 
distance or the photodissociation radius must be changed (see 
Sect. l8.2.Tl for a discussion of this). 

Figs[3]and|4]give the comparisons between the observed and 
the best-fit model CO lines for the carbon star LP And and the 
M-type star TX Cam, respectively. These results are discussed 
in Sects l8.1."Tl and [8T2~Tl respectively. The results for the remain- 
ing stars are shown in Appendix[C] The differences between the 
integrated line intensities in the best-fit model and the observed 
lines are presented in Table IC.l1 

7.2. Dust emission modelling 

7.2.1. Strategy for finding the best-fit model 

The observational constraints in the form of SEDs, covering the 
wavelength range ~1 -1300 /mi, are analyzed using ^-statistics. 
When combining the measured fluxes into an SED, the variabil- 
ity of the star must be taken into account. This is particularly 
important for fluxes measured at shorter wavelengths, where the 
amplitude of the variation is larger. To assess how strongly the 
fluxes are affected, the mean and deviation was calculated from 
several measurements performed at different epochs. As a result 
a 60% uncertainty was assumed for the J-band (1.24//m), 50% 
for the H-band (1.63/tm), 40% for the K-band (2.19/im), and 
30% for the L- and M-bands (3.79 and 4.64//m, respectively). 
At wavelengths longer than 5 /mi, a calibration uncertainty of 
20% was added to the total error budget and it dominates the er- 
ror (since there is limited information on the variability at these 
wavelengths). In some cases where there are many fluxes mea- 
sured at roughly the same wavelength, we introduce an average 
value to avoid a bias in the^ 2 calculation at certain wavelengths. 
(This is the reason why some of the fluxes presented in Table lB.ll 



are not included in the analysis.) The SEDs constructed from the 
photometric flux densities are thus averages, where the variabil- 
ity of the star at different wavelengths is reflected in the error 
bars (typically larger uncertainties at shorter wavelengths). The 
best-fit model is found from minimizing 

2 [ (-FmtxU _ fobs,,)) 1" 

where F,\ is the flux density and cr A the uncertainty in the mea- 
sured flux density at wavelength A, and the summation is done 
over all N independent observations. 

7.2.2. Results from dust emission modelling 

Table [5] lists the results from the dust modelling, the dust tem- 
perature at the inner radius T^fa), dust optical depth at lOjum 
T i0/jm> Afred f rom fitting the SED, and the number of observa- 
tional constraints N. ^ 2 -maps, at constant best-fit stellar ef- 
fective temperature, showing the level of consistency between 
the modelled SED and the photometric fluxes, when varying 
the dust temperature at the inner radius and the optical depth 
at 10 /mi, are displayed in Fig. [5] The SEDs from the dust- 
radiative-transfer model, with the observed photometric fluxes, 
are shown in Fig. [6] x 2 kA are generally close to 1 and the dust op- 
tical depth at 10 //m is rather well constrained for most sources 
(see Fig. |5l). 

7.2.3. IR spectra - a consistency check 

Low resolution spectra in the wavelength range 7.5-23 /mi cov- 
ering the silicate features at around 10 and 20 /mi, were mea- 
sured with the Low Resolution Spectrograph (LRS) as part of 
the IRAS mission. Later, the Short Wavelength Spectrometer 
(SWS) and the Long Wavelength Spectrometer (LWS) onboard 
ISO measured medium and high resolution spectra in the wave- 
length range 2.38-45.2/mi, and 43-197 /mi, respectively. SWS 
spectra also cover the silicate features and LWS spectra show 
the slope of the SED at longer wavelengths and is particularly 
sensitive to the cold dust. 

The strengths of the silicate features are v ery sensitive to th e 
dust optical depth [see for instance Fig. 1 in Ivan Loonl (120061) 1 
and in turn the dust mass-loss rate. As a consistency check, the 
IRAS-LRS spectra and SWS spectra are plotted for the M-type 
stars in the lowest panel of Fig. [6] The carbon stars all have the 
SiC feature at 11 //m (e.g.. ISpeck et alJll997h . However, since 
SiC grains compose a small fraction of the total dust content 
(~10%), and have a minor eff ect on the overall SED except for 
the 1 1 /mi feature (ISuhll2000h . we have not included SiC in the 
dust model and consequently the observed spectra for the carbon 
stars are not shown. LWS spectra are available for all stars but 
one, and they are also plotted in Fig. [6] 

A fit to the LRS and SWS spectra has not been attempted be- 
cause of the difficulty in including them in a proper best-fit anal- 
ysis. One complication is that the spectral shape and strength 
of the 9.7 /mi silicate feature can change with the pulsational 
phase i n stars with clear fe atures and definite bolometric vari- 
ability dMonnieret alj|1998l) . This is likely valid for all M-type 
stars in our sample. Also, since these spectra are reasonably well 
fitted (see Fig [6}, we expect the change in the resulting dust op- 
tical depth and inner radius, if the spectra were fitted in more 
detail (such as silicate features), to be small. A small change in 
the optical depth will result in an even smaller change in the es- 
timated mass-loss rate from the dynamical model (see Table|7]i. 
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Fig. 2. ^f 2 -maps showing the level of consistency between the velocity-integrated CO line intensity from the observations and the 
line model when varying the mass-loss rate and the h-parameter. The innermost contour corresponds to a 68% confidence level. The 
star indicates the best-fit model. The shapes of the line profiles were taken into account when choosing the best-fit model. 




Fig. 3. The observed CO J = 1 -> line from OSO and the 2 -> 1, 3 -> 2, 4 -> 3, and 6 -> 5 lines from JCMT (histogram) of the 
carbon star LP And, overlayed with the best-fit model (M=7 x 10~ 6 M yr -1 ) lines (solid line) from the radiative transfer modelling. 
The reduced^ 2 for this model is 0.25. See Sect. 18. 1 .Tl 
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Fig. 4. The observed CO J = 1 -> line from OSO and the 2 -» 1, 3 -> 2, and 4 -» 3 lines from JCMT (histogram) of the M-type 
star TX Cam, overlayed with the best-fit model (M=7 x 10~ 6 M Q yr -1 ) lines (solid line) from the radiative transfer modelling. The 
reduced^ 2 for this model is 3.0. See Sect. l8.2.T1 



For some stars (CW Leo, RW LMi, AFGL 3068, and IRC- 
10529) there is a slight discrepancy between the slope of the 
LWS spectrum and that of the best-fit model. To perfectly repro- 
duce the slope of the LWS spectra and the photometric fluxes 
at longer wavelengths, a varying mass-loss rate of these stars 
would probably have to be invoked. This h as already been sug- 
gested for CW Leo bv lGroenewegenld 19971) . who found that the 
dust density far from the star was underestimated by the best- 
fit r 2 model. Other possible causes could be a difference be- 
tween the optical properties of the model dust a nd the actu al 
dust surrounding the star [see e.g. Figs 1 and 6 in lSuhl d2000l) l. 



or the uncertainty in the long-wavelength dependence of the dust 
opacity. The fluxes beyond ~ 800 fim might also be affected by 
other processes than the dust emission, e.g., free-free emission 
and m olecular line emission dGroenewegenll 1 997t iDehaes et al.l 
120071) . Bearing in mind the uncertainties in the circumstellar 
model (Sect. [4}, and our aim to analyze all stars with the same 
assumptions we find that the LRS, SWS, and LWS data are well 
reproduced by our best-fit models. 
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Fig. 5. ^ 2 -maps showing the level of consistency between the photometric fluxes and the dust-continuum model when varying Tio^m 
and T&, at constant T e g. The innermost contour corresponds to a 68% confidence level. The star indicates the best-fit model chosen 
using the strategy described in Sect. I7.2.T1 All grids are terminated at T&{ri) = 1500 K since a dust temperature at the inner radius 
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Fig. 6. Observed photometric fluxes together with ISO-LWS and ISO-S WS spectra overlayed by the best-fit model SEDs derived 
from the dust radiative transfer model. The model is represented by the solid line in all panels. The dashed line shows the model 
SED when the interstellar extinction is not taken into account. The lower panels show the wavelength interval 5-35 /im in more 
detail for the M-type stars. The best-fit models are here shown together with the observed IRAS-LRS spectra (dotted line) and the 
ISO-SWS spectra (dot-dashed line) when available. The ISO-LWS spectra are plotted for all stars but GX Mon. 

7.2.4. Mass-loss-rates from dust continuum and dynamical Eqns. © and © corresponds to the observed expansion veloc- 
model ity within 0.1%. 



In the dynamical model, the density structure is calculated from 
the conservation of mass and the model is iterated until this 
gives an optical depth that agrees with the dust continuum opti- 
cal depth obtained from the SED modelling (within one percent) Table[6]lists mass-loss rates, i'd r , T*, and (Qrp) ' f° r the best-fit 
and the resulting terminal gas expansion velocity from solving dynamical model. 
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Table 5. Results from the dust emission modelling. 



Source 




7d(n) 


no 


*r 2 ed 


N 




[K] 


[K] 






LP And 


2400 


1100 


0.55 


1.6 


11 


CW Leo 


2000 


1500 


1.2 


2.4 


12 


RWLMi 


2400 


1200 


0.6 


0.9 


1 1 


V384 Per 


1800 


1000 


0.2 


0.9 


10 


AFGL 3068 


2400 


1500 


4.0 


5.3 


11 


TX Cam 


2200 


800 


0.4 


0.5 


9 


GX Mon 


1800 


500 


0.3 


4.9 


9 


WXPsc 


1800 


1000 


3.0 


1.1 


12 


IKTau 


2400 


1200 


1.2 


0.9 


12 


IRC- 10529 


1800 


800 


3.0 


3.1 


10 



lines have resulted in mass-loss rates close to our results , 
1.2 x lO^Moyr 1 (C O/FW.l x 1(T 3 ) jTevssier et all [2006b. 
and 1.5 x 10" 5 M yr 1 dSchoier & Olofssonll200ll). b oth assum- 
ing a distance of 630 pc. iMen'shchikov et alj d2006l) estimated 
the mass-loss rate to be 1.9 x 10~ 5 M Q yr~' from a dust radiative 
transfer model, assuming a distance of 740 pc and a dust-to-gas 
mass ratio of 0.0039. The corresponding mass-loss rate at the 
distance used and the dust-to-gas ratio obtained in our model, 
would be 1 . 1 x 10~ 5 M Q yr , in good agreement with our result. 

iNeri et al.l d!998l) found a slight asy mmetry in the CSE 
from m aps of the CO(/ = 1 — > 0) emission. Mauron & Hugginsl 
(120061) . using optical images in scattered light, suggested that 
the envelope around LP And may have an interesting structure, 
but concluded that the signal-to-noise in their data is too low to 
make any firm statement. 



Table 6. Mass-loss rates and other quantities from the dynamical 
model. 



Source 


M 


Vdr 






M d y n /Mc0 


[io- 


5 Moyr 


[kms- 1 ] 


[io- 3 ] 






LP And 


17 


3.1 


1 


0.09 


2.5 


CW Leo 


21 


1.9 


2 


0.05 


1.1 


RWLMi 


17 


1.3 


2 


0.08 


3.3 


V384 Per 


5 


6.5 


2 


0.11 


1.6 


AFGL 3068 


62 


0.7 


3 


0.02 


3.1 


TXCam 


3 


10.6 


4 


0.11 


0.5 


GX Mon 


2 


12.6 


10 


0.09 


0.1 


WXPsc 


18 


4.3 


3 


0.06 


0.5 


IKTau 


7 


7.3 


2 


0.10 


0.7 


IRC-10529 


26 


2.9 


2 


0.06 


1.6 



8. Individual objects 

All sources in the sample are well-studied both at IR and radio 
wavelengths. As established already in Sects l7.ll and l7.2l we are 
able to reasonably successfully model the CO line emission, the 
integrated intensities and the line shapes, and the SEDs for all 
10 stars using a constant mass-loss rate. 

In the following discussion it should be kept in mind that the 
different methods measure the mass-loss rate at different epochs 
and over different time intervals. The size of the CO line emit- 
ting region decreases with the /-number and hence the different 
lines measure the mass-loss rate at slightly different epochs. The 
same is true for the SEDs where the emission at different wave- 
lenghts (apart from the photospheric contribution) comes from 
different parts of the CSE. The estimated mass-loss rate should 
be considered as an average and it is a complicated synthesis of 
the different epochs and time scales. 

8. 1 . Carbon stars 
8.1.1. LP And 

The derived mass-loss rate from the CO modelling is 
7 X lO^Moyr -1 , while the dynamical model gives a mass-loss 
rate of 1.7xlO- 5 M yr a factor of 2.5 larger. The model 
lines agree very well with all five observed lines, ;^ ed =0.25 
(Fig. [3]). Previous mass-loss-rate estimates from CO radio 



8.1.2. CW Leo (IRC+10216) 

Being the closest carbon star, CW Leo is certainly the most well- 
studied. We estimate the mass-loss rate to be 2 x lO^Moyr -1 
from the CO line modelling and 2.1 x lO^Moyr 1 from 
the dynamical model. Previous es timates of its mas s-loss 
rate from CO is 1. 2 x 10' 5 Mm vr' 1 (fTeyssier et al.l 120061) , and 
1.5 x lO^Mgyr 1 dSchoier & Olofssonl 1200 lb. bo th assuming 
a distance of 120pc. Groenewegen et alT i 19981) did a very 
thorough investigation of the CO radio line emission toward 
CW Leo and estimated the present day mass-loss rate to 
(1.5+0.3) X 1O- 5 M yr -1 and the distance to lie between 1 10 and 
135 pc. 

Several reports on structures in the circumstellar environ- 
ment around CW Leo c a n be found in th e literature (e.g., 
iMauron & Hugginsl 119991: iFong et alj 120061) . Optical images 
show arc-like structures on arcsecond scales in the CSE, but 
whether this is due to chan ges in the mass-los s rate over time, 
is still under debate [e.g.. iLeao et alj d2006l) and references 
therein]. 



8.1.3. V384 Per 

The best-fit CO model gives a mass-loss rate of 
3 x 10~ 6 M G yr -1 . However, as can be seen in Fig. IC.3I the 
model y = 1 — > 0, 2 — > 1 , and 4 — > 3 lines are all 15-20% weaker 
than the observed lines, while the model J = 3— >2 line is 65% 
stronger than the observed line. This results in a rather high 
;^ 2 ed =2.8. The most likely explanation is a calibration error 
in the / = 3 — » 2 line. The mass-loss rate from the dynamical 
model is higher, M = 5 X lO^Mgyr" 1 . ISchoier & Olofssonl 
( 1200 lb estimated M = 3.5 x lO" 6 M yr _1 at 560 pc from CO 
line modelling. 



8.1.4. AFGL 3068 (LL Peg) 

Our model of the CO line emission reproduces all five lines 
for a constant mass-loss rate of 2 X 10~ 5 M yr -1 with a 
xLa=0.7. The dynamical model gives M = 6.2x 10~ 5 M o yr _1 , 
slightly more than a factor of three higher. Earlier estimates 
of the mass-loss rate fro m CO line modellin g have resulted 
in M = 6.0 x IO 5 M^vr" 1 JTevssier et _ al.| [2006b a t lOOOpc, and 
2.0 x l O^Meyr 1 dWoods et al.ll2003b at 820 pc. iTevssier et al.l 
d2006b deduced a much lower temperature in the outer envelope 
than we calculate in our model and this is the reason for the dis- 
crepancy between the estimates. 
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There is no doubt that the CSE around this source is asym- 
metric and have density struct ures that cannot be reconci led with 
a spherically symmetric CSE dMauron & Hugginsl2006l) . Radial 
CO (J = 1 — > 0) visibility profiles show a compact inner enve- 
lope surrounded by a large outer shell. The outer shell is strongly 
bipolar and there are even indications that it might be detached 
dNeriet al.lll998l) . 

8.1.5. RW LMi (CIT 6) 

We find the CO line data to be consistent with a mass-loss 
rate of 5 x 10~ 6 M yr _1 . From the dynamical m odel the mass- 
loss ra te is estimated to be 17 x 1(T 6 M yr _1 . iTevssier et alj 
(120061) found it necessary to infer a change in the mass- 
loss rate about 250 years ago to fit the CO intensity mea- 
sured at offset positions; the mass-loss rate has decreased from 
7.5 x 10 _6 M Q yr _1 to 5.0 x 1O~ 6 M yr' 1 . The two component 
structure was already noticed by N eri et al.l d!998l) when first 



inspecting the interferometric data. Schoi er & Olofssonl d200ll) 
estimated M = 6.0 x 10~ 6 M G yr~' f r om CO line modelling using 



the same distance as ITevssier et al.l d2006l) . 440 pc. 

8.2. M-type stars 
8.2.1. TX Cam 

The best-fit CO mo del gives M - 7 x 10~ 6 M yr~' [the same 
rate was derived by iGonzalez Delgado et al.l (|2003|) 1. but the 
model does not reproduce the shape of the observed line profiles 
(Fig|4]i. The observed J — 1 — > and 2 — > 1 lines are fiat-topped, 
which suggests unresolved, optically thin emission, while the 
model gives double-peaked lines, which suggests resolved, opti- 
cally thin emission. The observed ./ = 3 — > 2 and 4 — > 3 lines are 
parabolic and thus optically thick, while the model produces an 
optically thin, resolved J — 3 — > 2 line and a flat-topped 4 — > 3 
line. The simplest explanation would be that the distance is un- 
derestimated. A larger distance requires a higher mass-loss rate 
to reproduce the line intensities and hence gives line profiles 
that better reproduce the observed ones. An increased distance 
also decreases the angular size, but this is partly counteracted by 
the fact that the CO photodissociation radius increases with the 
mass-loss rate. Another explanation could be that the size of the 
CO envelope in this object is smaller than the phot odissociation 
radius given by the model of iMamon et al.l d!988l) . Decreasing 
the size of the CO envelope solves the problem with the model 
producing resolved lines. Also, a smaller radius will require a 
higher mass-loss rate to reproduce the observed intensities. This 
also increases the optical depth. However, the dynamical model 
estimate of the mass-loss rate is 3 x lO _6 M yr _1 , a little more 
than a factor of 2 less than the CO modelling. 



locity. [Setting T&(r\) = 1000 K, requires a higher optical depth 
to fit the model to the observed fluxes (tio=0.85), and results 
in M — 6.3 x 10~ 6 M yr _1 from the dynamical model, somewhat 
more in line with the CO results. The measured expansion ve- 
locity would be reached with *F as 8.3 x 10~ 3 .] A possible ex- 
planation for the low temperature at the inner radius is that the 
star has undergone a rather drastic decrease in the mass-loss rate 
recently. The dust envelope from an earlier epoch would have 
moved outward and cooled to the temperature obtained from 
the modelling, and the photometric fluxes would thus probe this 
'older' shell. In such a case the dynamical model would produce 
an erroneous result. 



8.2.3. WX Psc 

The / =1— > 0, 2— > 1, 3 — > 2, and 4 — > 3 lines are well-fitted by a 
constant mass-loss rate of M = 4x 10~ 5 M yr _1 . We are, how- 
ever, not able to fit the ./ = 6 — > 5 line using this model. The 
problem of not being able to model higher-frequency CO lines 
from this sourc e using a cons t ant ma s s-loss rate has been p revi- 
ously noted by Kemper et al.l d2003l) . ITevssier et alj d2006l) and 
iDecin et ail d2007l) . The difficulty of getting well-calibrated data 
for these higher frequency lines should be kept in mind when 
drawing conclusions based on this kind of data. As an example, 
the integrated intensity in the 6 — » 5 li ne in this work is a f ac- 
tor of two stronger than that reported bv lKemper et al.l d2003l) . A 
model fitting only our 6 — > 5 line results in M = 4 x 10~ 6 M yr~' . 
IDecin et alj d2007l) introduce a rather complicated mass-loss his- 
tory, and estimate the mass-loss rate to be 6 x 10~ 6 M yr~' in the 
region where the 6 — > 5 line is mainly coming from. Higher tran- 



sition CO lines are within the range of archived ISO LWS spectra 
(43 - 197 fim). The CO J = 16 -> 15 and 7=18^17 lines were 
modelled using our best-fit model. The resulting peak intensities 
[4.6 x lO^Wcm-^km/s)- 1 and 5.0x 10~ 22 WctrT 2 (km/s)- 1 , 
respectively] would result in a S/N a 5 [the resolution-corrected 
lo- noise level is 8.9 x 10~ 23 W cirT 2 (km/s)" 1 ]. None of the CO 
lines within the range of the ISO LWS were detected, and this is 
in line with the weakness of the observed J — 6 — > 5 line. Thus, 
the results reported by us and previous authors suggest a recent 
decrease in the mass-loss rate in this source. From the dynamical 
model we find M= 1 . 8 x 1 0~ 5 M yr~ 1 . 



8.2.4. IK Tau 

Modelling IK Tau with a constant mass-loss rate of 
1 x 10~ 5 M yr _1 gives a good fit to the four available CO lines 
(Fig.lC.7[>. Th i s rate is a factor of two higher than that found by 



ITevssier et alj d2006l) at a somewhat smaller distance of 250 pc 
The dynamical model gives M=7x lO^Moyr 1 



8.2.2. GX Mon 

The dynamical model gives a mass-loss rate estimate 
(2x lO _6 M yr _1 ) one order of magnitude lower than the CO 
model (2xlO" 5 M yr 1 ). The best-fit SED model found ac- 
cording to the strategy described above gives a dust temper- 
ature at the inner radius of 500 K. Fig. [5] shows the ^ 2 -maps 
from the SED-fitting and in the case of GX Mon, the model is 
not very well constrained. The major dust condensation phase 
in M-type stars is thought to occ ur at temperatures between 
1200K and 800 K (I Whittetl [2003b . With the low dust optical 
depth derived for GX Mon a high dust-to-gas ratio Q¥ = 0.01) 
is needed to accelerate the gas to the measured expansion ve- 



8.2.5. IRC-10529 

When modelling IRC-10529, we encounter a similar problem 
as for V384 Per. We find that the best-fit CO model with 
M = 3 x 10~ 5 M yr~ 1 has a * 2 ed =3.8. The J = 1 -» and 3 -> 2 
lines are fitted within the observational uncertainties, but the 
model 2 — > 1 line is 45% too weak, and the 4 — > 3 line is 40% 
too strong (Fig IC.8I ). We do not find it likely that this is due 
to modulations in the mass-loss rate, since the emitting re- 
gions for the 1 — > and 2 — > 1 lines overlap to a large extent. 
It is more likely due to observational problems. The dynamical 
model gives M = 2.6 x 10 5 M yr _1 , in good agreement with the 
estimate from the CO model. 
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Table 7. The resulting change in mass-loss rate (M) and dust-to- 
gas ratio (W) when varying the input parameters in the dynamical 
model ofRWLMi. 



Parameter 


Change 


M 


*P 


L 


+50% 


+ 50% 


- 35% 




-50% 


- 50% 


+ 100% 


M 


+50% 


- 25% 


+ 30% 




-50% 


+ 30% 


- 30% 


0"d 


+50% 


+ 50% 


- 35% 




-50% 


- 50% 


+125% 


»'d 


+50% 


- 35% 


+ 60% 




-50% 


+ 105% 


- 55% 


TlO 


+50% 


+ 35% 


+ 10% 




-50% 


- 40% 


- 10% 


Qa 


+50% 


+ 50% 


- 35% 




-50% 


- 50% 


+ 100% 



9. Discussion 

9.1. Dependence on parameters 

The model CO line intensities dependence on the /^-parameter 
and M is displayed in Fig. [2] A further discussion on how the 
resulting li ne intensities depend on vari ous input parameters can 
be found in lSchoier & Olofssonl(l200ll) . 

To evaluate the sensitivity of the results of the dynamical 
model to the different derived or assumed input parameters, we 
varied the parameters in Eqns dUl and (© by + 50%. Shown in 
Table[7jis the needed change of the dust-to-gas ratio, VP, and the 
mass-loss rate, M, to reach the observed expansion velocity and 
derived optical depth for RW LMi. As can be seen, the dynam- 
ical model results are strongly dependent on the assumed dust 
parameters, which are not well known. 

A trend with mass-loss rate can be seen in the dynam- 
ical model. For TX Cam, an intermediate-mass-loss-rate ob- 
ject (M = 3 x 10~ 6 M yr~' from the dynamical model), an in- 
crease in the mass-loss rate by 50% leads to an increased expan- 
sion velocity by 6%. For RW LMi, a high-mass-loss-rate object 
(M = 1.7 x 10~ 5 M yr~') the same increase in the mass-loss rate 
will lead to a decrease in the expansion velocity by 15%. This is 
due to the fact that a higher mass-loss rate leads to a higher dust 
optical depth resulting in a radiation field that decreases more 
rapidly with radius and with it the radiation pressure efficiency. 
Nevertheless, this does not lead to a problem finding the correct 
mass-loss rate. In some cases, two mass-loss rates give the same 
expansion velocity, but there is only one solution that produces 
the correct dust optical depth. 

We have decided not to display any error bars on the esti- 
mates from the dynamical model in Fig. UJ The reason for this 
is that the resulting (formal) error bars given by the dust radia- 
tive transfer model (variations in inner radius and optical depth, 
Fig. [5j are small, and ignores the substantial uncertainties in the 
adopted dust parameters. 

9.2. Comparing the two methods 

It is very difficult to assess which method is the most reliable for 
estimating mass-loss rates, since both are dependent on a num- 
ber of assumptions, which are more or less well-founded. For 
instance, the dynamical model assumes an entirely dust driven 



wind, while recent results suggest that there might b e other fac- 
tors at play dWoitkdl2006tlHofner & Andersenll2007l) . 

In all cases but GX Mon, for which we have good reasons 
to believe that the dynamical model gives an erroneous result 
(see Sect. 18.2.21 for an explanation), the mass-loss rates esti- 
mated from the dynamical model and the CO line modelling 
agree within a factor of ~ 3. A comparison is shown in Fig. [7j 
where the full drawn line marks the one-to-one relation, while 
the dashed lines subtend the range within a factor of 3. [Error 
bars represents the range within the innermost contour in the^f 2 - 
maps from the CO model. Here the shape of the line profiles have 
also been taken into account.] Since the two methods give inde- 
pendent estimates of the mass-loss rate, this result can be used 
as a guideline to the uncertainty in the mass-loss-rate es t imates 
within the adopted circumstellar model. iHeras & Honvl d2005[) 
found in their comparison of mass-loss rates derived from ISO- 
SWS and CO mass-loss rates taken from the literature for low- 
mass-loss-rate M-type stars that the estimated mass-loss rates 
can differ up to a factor of ten. 

In principle, it is possible to check whether the dust results 
obtained from the dynamical model (v<jr> *P, (Grp) f ) are consis- 
tent with those used in the CO radiative transfer modelling (h, 
(Qrp) ). However, we do not expect to get a perfect agreement 
since in the CO radiative transfer modelling the /z-parameter is 
adjusted until the radial kinetic temperature profile results in rea- 
sonable relative strengths of the CO lines. Hence, it is affected by 
possible missing heating and cooling terms in the description of 
the thermodynamics of the circumstellar gas. However, as out- 
lined in Sect.[5j this does not affect the mass-loss-rate estimates. 

For the carbon stars in our sampl e the dust-to-gas ratio, is 
estimated to be 0.001 < *F < 0.003. iGroenewegen et ail d 19981) 
derived dust-to-gas ratios for a sample of 36 carbon stars and 
found a rather constant value of 0.0025 up to P « 500 where 
a strong increase with period is found. This agrees rather well 
with our results. However for AFGL 3068 they find *F * 0.02, in 
strong excess of our result. For the M-type stars we find some- 
what large r dust-to-gas rat i os, 0.0 02 < *F < 0.004 (GX Mon 
excluded). lJusttanont et alj d 19941) modelled three thick M-type 
CSEs and found 0.003 < HP < 0.006. 

The CO mass-loss-rate estimate is lower than the estimate 
from the dynamical model for all five carbon stars, see Fig. |7J 
The opposite can be seen for the M-type stars, where the CO 
model gives a mass-loss-rate estimate higher than the dynamical 
model for three out of four M-type stars (GX Mon excluded). 
Keeping in mind that the sample presented here is too small to 
make a general statement, we suggest that a possible explana- 
tion to this is that the adopted CO/H2-ratio is too low, on aver- 
age, for the M-type stars and too high for th e carbon stars. Thi s 
is also in line with recent chemical models dCherchneffl 12006 ) . 
where a non-equilibrium chemical model gives a CO-abundance 
of 6 x 10" 4 in M-type stars (C/O=0.75) at 5 stellar radii. For the 
carbon stars (C/0=1.1), the resulting CO-abundance at the same 
distance from the star is found to be 9 x 10 4 . Adopting these 
abundances would require a higher mass-loss rate for the carbon 
stars and a lower mass-loss rate for the M-type stars to reproduce 
the observed CO line intensities. 

9.3. Comments on some results by Kemper et al. (2003) 

iKemper et al.l d2003l) made non-LTE radiative transfer modelling 
of CO / = 2 — > 1 to 7 = 6— >5 emission toward a sample of AGB 
stars and red supergiants. In a line ratio diagram they compared 
modelled line intensity ratios, for standard CSEs with different 
mass-loss rates, to those found from observations (Fig. 4 in their 



12 



S. Ramstedt et al.: On the reliability of mass-loss-rate estimates for AGB stars 




Fig. 7. Comparing the mass-loss-rate results from the CO model 
and the dynamical model. Dots mark the M-type stars, and 
squares the carbon stars. The full drawn line represents a one- 
to-one relation and the dashed lines subtend the range within a 
factor of 3. 




1 2 
I(3-2)/I(2-l) 

Fig. 8. Observed line intensity ratios for the sample sources 
(squares) compared with predictions from the standard model 
for various mass loss rates (stars). 



paper). Kemper et al. found that they could not reproduce the ob- 
served line intensity ratios. On the contrary, the standard model 
described above (Sect. |5j, generally gives a very good fit to the 
observed line intensities as indicated by a Xied of the order of 
unity for the majority of the sample sources. This is also illus- 
trated in Fig. [8] where line intensity ratios calculated from our 
standard model using v e = 15 km s and h = 0.5 for a range of 
mass-loss rates are compared to the observations. The model 
CSE was placed at a distance of 1 kpc. The interconnected stars 
are the results from the standard model and the squares show the 
ob served line rati os. 

iKemper et alj d2003l) also found that the magnitude of the 
stochastic velocity strongly affects the line intensities, since 
it determines the size of the region where the line is formed 
[v st0 = (vf urb + v^ erm ) - 5 ~ v,urb, since v, herm is the local ther- 



mal broadening and v t h e rm VtmbL We have performed the 
same test as in their Fig. 10, changing the turbulent velocity 
from 0.05 to 2 km s _1 , to see how this affects the intensity of the 
CO(7 = 3 — > 2) line. We used the same standard model as above 
with a mass-loss rate of M-l x 10~ 5 M yr _1 . IKemper et al.l 
(120031) saw a decrease in intensity from around r m b=2.4K at 
v st o=0.05kms _1 to r m b=0.2K at v sto =2.0kms~ 1 . We, however, 
see the opposite trend and a much weaker effect; r m b=1.05K 
at v tU rb=0.05kms _1 to r m b=1.28K at v tur b=2.0kms _1 . It is dif- 
ficult to outline exactly how the line intensity is affected by 
the turbulent velocity, as it also affects the excitation of the 
molecules. However, to a first approximation, with an increased 
turbulent width the number of molecules per unit velocity inter- 
val decreases. On the other hand, as already mentioned above, 
it also increases the size of the interacting region. Under certain 
circumstances (spherically symmetric CSE, constant excitation 
temperature, v tu ,-b <K v e ) these two effects cancel out and the line 
intensity becomes independent of the turbu lent velocity (see for 
example lMorrill975l;loTofsson et al.lll982l) . 

In concl usion, we are not a ble to reproduce the results 
reported by IKemper et al.l d2003l) using our model. We con- 
firmed our result usin g two additional codes; an acc eler- 
ated Monte Carlo code (H ogerheiide & van derjak 2000 ) and 
an A ccelerated Lambda Iteration (ALI) code dMaercker et alj 
120081) . The results from the excitation analysis were then com- 
bined with different r ay tracers [the one used in t his pa per and 
the code described in lHogerheijde & van der Tiki (I2000M to get 
the resulting line profiles. All results were in good agreement 
with ours. 



10. Conclusions 

We have presented mass-loss-rate estimates, derived using two 
independent methods, for a sample of 10 intermediate- to high- 
mass-loss-rate, well-studied AGB stars, and arrive at the follow- 
ing conclusions: 

- The mass-loss rate estimates obtained using the two different 
methods agree within a factor of ~ 3, suggesting that this is 
(at least) the uncertainty in present mass-loss-rate estimates, 
e.g., using CO radio line emission. The average of the ratio 
between the two mass-loss-rate estimates is (Mdyn/^co) = 
2.3 + 1.0 for the carbon stars, and (M dyn /M C o) = 0.8 + 0.5 
for the M-type stars (excluding GX Mon). 

- In most cases, it is not possible to set an upper limit to the 
mass-loss rate for the high-mass-loss-rate objects by mod- 
elling only the integrated CO line intensities. The analysis 
can be improved by taking the shape of the line profile into 
account. 

- We are able to successfully model the CO(7 = 1 — > to 
4 — > 3) data of all stars using a constant mass-loss rate. 
However, in most cases there are various indications that our 
'simple' circumstellar model does not fully reproduce the 
data. Effects of time-variable mass loss, and non-isotropic 
mass loss may be present. 

- We do not find any trends with mass-loss rate when compar- 
ing the two models, but we find a possible dependence on 
chemical type. Generally, the CO model estimates a lower 
mass-loss rate for the carbon stars than the dynamical model, 
while the opposite behaviour can be observed for the M- 
type stars. We conclude that a CO/H2-ratio increased and 
decreased by about a factor of 1 .5 for the M- and C-stars, 
respectively, would give a very good agreement between the 
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methods, but caution that this is a tentative conclusion con- 
sidering that the sample is small, and there may be other sys- 
tematic errors. 
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Appendix A: A new mass loss rate formula 

iKnapp&Morrisld 19851) derived a theoretical mass-loss-rate for- 
mula in terms of the CO(/ = 1 — > 0) line intensity, valid for op- 
tically thick CO emission. The formula relates the mass-loss 
rate to a number of easily determined observables; T m b for the 
CO(l-O) line, and v e . The distance D and CO abundance fco 
also enter the formula. It was derived u sing a fixed CO enve- 
lope size. This formula was used in e. g. lOlofsson et al.l d 1993b 
to estimate mass-loss rates for a sample of carbon stars. Later, 
ISchoier & Olofssonl d200ll) could compare their results on the 
sample, derived using the same radiative transfer model as in 
this work, to those estimated using the formula. They found 
that the formula underestimated the mass-loss rates compared to 
what was found from the radiative transfer analysis, on average 
by a factor of four, but the discrepancy was found to increase 
even more for lower mass-loss rate objects. The reason is that 
the emission from the low-mass-loss-rate objects is not optically 
thick and it is more sensitive to the size of the CO envelope. 
However, the usefulness of simple formulae of this sort can not 
be disputed. 

In order to improve the results by iKnapp & Morris! d 19851) . 
we set up a grid with 60 model stars, whose properties 
were varied in M (1, 3, 10, 30, and 100 x 10~ 7 M yr" 1 ), v e 
(5, 10, 1 5 and 20kms 1 ), and / C p (1, 3, and lOxlO" 4 ). 
Following ISchoier & Olofssonl feOOll) . the value of h was as- 
sumed to be 0.2 for M in the range up to 3 x 10~ 7 M Q yr -1 , 
0.5 for 1 - 3 x 10 _6 M o yr _1 , and 1.5 for the mass-loss rate 
1 x 10~ 5 M Q yr -1 . For each star, the velocity integrated inten- 
sities in CO (/=l->0, / = 2-»l, / = 3->2, / = 4-»3) were 
calculated for a 20 m class telescope (the corresponding beam 
sizes were 33", 16.5", 11", and 8.5"). The distance to the stars 
were set to 1 kpc. The output from the modeling was fitted to 
Eqn. dA.lt and the errors were minimized using a Levenberg- 
Marquardt method (Press et al. 1992), 



CO 



(A.l) 



Table I A. 1 1 summarizes the results. For low upper-/ levels, 
the mass-loss rate is about equally (and moderately) dependent 
on changes of /, v e , and fco. At high upper-/ levels, the mass- 
loss rate is almost linearly dependent on the velocity-integrated 
intensity (/), but relatively insensitive to changes in v e , and /co- 
in general, the parameter most difficult to fit was the expo- 
nent for the expansion velocity (v e ), while the exponent for the 
velocity-integrated intensity (/) is well determined. Eqn. dA.ll i 
only applies to unresolved emission. Fig. IA. II shows a compar- 
ison between mass-loss rates estimated using Eqn. dA.lt on the 
y-axis and mass-loss rates estimated using the Monte-Carlo ra- 
diative transfer model on the x-axis for three samples of C- , 
M-, and S-stars dSchoier & Olofssonl200lHOlofsson et al.l2002t 
iGonzalez Delgado et all 120031: iRamsted t et al. 20061). f co is as- 
sumed to be 1 X 10~ 3 for the C-stars, 2 x 10~ 4 for the M-stars, 
and 6 x 10~ 4 for the S-stars. 




10 8 10 7 10 6 10 5 10 4 10 8 10 7 10 8 10 5 io ■ 




Fig. A.l. Mass-loss rates derived using Eqn. lA.ll and the param- 
eters in Table lA.ll versus the mass-loss rates estimated using the 
Monte Carlo model for the same stars. Red triangles represent 
carbon stars, blue squares, M-type stars, and green dots repre- 
sent S-type stars. 



Table A.l. Values for the parameters in Eq. IA.ll sj is calibrated 
such that Eq. lA. II gives the mass-loss rate in M yr -1 . Values are 
valid in the ranges 10~ 7 <M< lO~ 5 M yr~', 5 < v e <20kms~', 
10~ 4 < fm < 10~ 3 . 
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bj (V e ) 


cj (fco) 


/ = 


l->0 
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t0.03) 
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0.59(±0.07) 
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Appendix B: Photometric fluxes and dust radiative transfer 

The SEDs used as constraints in the dust emission modelling are constructed from JHKLM-band, sub-millimetre, and millimetre 
photometric flux densities from the literature, together with IRAS data. All flux densities used are listed in Table IB. 11 Values in 
italics were omitted from the analysis to avoid a bias in the ^-analysis at certain wavelengths. Values in italics marked with a 
dagger were omitted since the reported errors were > 50%. 



Table B.l. Photometric flux densities in 10 26 W m~ 2 Hz" 1 . 



wavelength 
[jum] 


LP And 


CW Leo 


V384 Pel- 


AFGL 3068 


RWLMi 


TX Cam 


GXMon 


WX Psc 


IK Tau 


IRC - 10529 


1.24 


1.6 a 


4.0 d 


ls' 5 




5.8 h 


280" 


49 s 


3.2' 


250' 


2.4' 


1.63 


IP 


58 d 


53 c 




39 h 


590 a 


120) 


32' 


730' 


20' 


2.19 


55 a 


520 d 


170 e 




180 h 


780 a 


200* 


130 1 


1200 1 


78' 


3.79 


450 a 


10000 d 


520 e 


9.68 


970 b 


750 a 


580* 


620' 


1900' 


290' 


4.64 




19000 d 




55 s 






410) 


800' 


1800 1 


340' 


8.55 








290 k 








&50 k 






8.57 
















700 k 






9.56 








360 k 








1100 k 






11.21 








40& 


3900 k 






1000 k 






12 


960 b 


48000 b 


540 b 


1700 b 


3300 b 


1600 b 


1600 b 


1200 b 


4600 b 


1300 b 


12.45 








47& 








360 k 






25 


470 b 


23000 b 


200 b 


780 b 


1200 b 


640 b 


360 b 


970 b 


2400 b 


1100 b 


60 


110 b 


5700 b 


40 b 


250 b 


270 b 


130 b 


110 b 


220 b 


330 b 


220 b 


100 


36 b 


920 b 


12 b 


74 b 


86 b 


39 b 


40" 


72 b 


100 b 


64" 


400 
















3.0 m 






450 


0.35 f ' 


13" 




2.2° 


2.4" 








1.2 H 




600 


0.41 f 




0Ji ft 


0.67 f 














800 


0.25 f 


6.7 1 


0.21 f 


0.33 1 


0.78 1 








0.25 1 




850 




12" 




0.59° 


0.97° 












1100 


0.081 f 


1.9' 


0.045 f 


0.17 f 


0.16 1 








0.082 1 




1200 
















0.24 c 




0.10 c 


1300 




1.5" 






0.29" 






0.059 c 


0.037" 





References: 
a |Bagnuldjl996h 

b Beichman etTTl dl988h . IRAS Point Source Catalogue 

c Dehaes et al] j2007h 

d Epchtein et al] Jl985h 

e GezarietaLH 19871) 

f Groenewegen et al] l l 19931) 

g Jones etal] i 19901) 

Kerschbaum & Hronl ( ll994h 
' iLeBertrelll 19931) 
j lLepineetal]j 19951) 
k lMarengo et al] (119991) 
' iMarshalletaiTj^^ 
m ISopkaetal.Ul985l) 
" IWalmslevetal] ( ll99lh 
° JCMT archival data 
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Appendix C: Results from the CO line radiative transfer model 

Figs IC.UC.8l show the best-fit models from the CO radiative transfer modelling overlayed the new observational data for all stars in 
the sample (except LP And and TX Cam shown in Figs[3]and|4j respectively). The results are discussed in Sect. [8] Table IC.fl gives 
the difference between in integrated intensity from the model and the observation, for each line. 
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Fig. C.l. The observed CO J = 1 -> line from OSO, 2 -> 1, 3^2, and 4 -» 3 lines from JCMT, of the carbon star CW Leo, 
overlayed with the best-fit model (M=2 x 10~ 5 M yr~') from the CO radiative transfer modelling. The reduced^ 2 for this model is 
xl d =0-6- See Sect. [8X2] 





- RW LMi 


C0(l-0)- 




OSO 






1 

EH 








O 









-50 

v [km s ] 




_ RW LMi 


C0(3-2) H 


_ JCMT 




50 





V 


[km s ] 


3- 


-> 2, and 4 -> 



~ RW LMi 


C0(4-3)^ 


JCMT /fr\ 











v [km s 



overlayed with the best-fit model (M=5 x 10 M yr ) from the CO radiative transfer modelling. The reduced^- for this model is 



^red" 


=2.1. See Sect.|8.1.5| 
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Fig. C.3. The observed CO / = 1 — > line from OSO, 2 -> 1, 3^2, and 4 -> 3 lines from JCMT, of the carbon star V384 Per, 
overlayed with the best-fit model (M=3 x 1CT 6 M Q yr~ [ ) from the CO radiative transfer modelling. The reduced^- 2 for this model is 
*L=3.7. See Sect. [TO] 
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Fig. C.4. The observed CO J = 1 line from OSO, 2 -» 1, 3 -» 2, 4 -> 3, and 6 -> 5 lines from JCMT, of the carbon star AFGL 
3068, overlayed with the best-fit model (M-2 x 10~ 5 M yr -1 ) from the CO radiative transfer modelling. The reduced^ 2 for this 
model is^ 2 ed =0.7. See Sect. l8X4l 
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Fig. C.5. The observed CO J = 1 -> line from OSO, 2 -» 1, 3 -> 2, and 4 -> 3 lines from JCMT, of the M-type star GX Mon, 
overlayed with the best-fit model (M-2 x 10~ 5 M yr -1 ) from the CO radiative transfer modelling. The reduced^ 2 for this model is 
;tf ed =0.04. See Sect. 
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Fig. C.6. The observed CO J = 1 -> line from OSO, 2 -> 1, 3 -» 2, 4 -» 3, and 6 -> 5 lines from JCMT, of the M-type star WX Psc, 
overlayed with the best-fit model (M=4 x 10 5 M yr _1 ) from the CO radiative transfer modelling. The reduced^ 2 for this model is 
^- 2 d =1.0, when the 6 — >5 line is omitted. See Sect. 18.2.31 
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Fig. C.7. The observed CO J = 1 -> line from OSO, 2 -> 1 , 3 -> 2, and 4 -> 3 lines from JCMT, of the M-type star IK Tau, overlayed 
with the best-fit model (M=l x 10 _5 M o yr _1 ) from the CO radiative transfer modelling. The reduced^- 2 for this model is;f 2 ed =0.5. 
See Sect. [8^41 




Fig. C.8. The observed CO / = 1 — > line from OSO, 2 -> 1, 3 -> 2, and 4 -> 3 lines from JCMT, of the M-type star IRC- 10529, 
overlayed with the best-fit model (M=1.6 x 10~ 5 M o yr -1 ) from the CO radiative transfer modelling. The reduced^- 2 for this model 
is^ 2 ed =4.6. See Sect. l8~2~5l 
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Table C.l. The difference between the integrated CO(/ — > / - 1 ) intensities from the best-fit model and the observations, in percent; 

(^mod ^obs)/^obs- 



Source 


OSO 




JCMT 






1 


2-» 1 


3 — > ? 


4 — > 3 


6^5 


LP And 


3.2 


5.4 


15.8 


4.2 


-0.6 


CW Leo 


4.8 


-7.3 


-15.6 


-11.7 




RWLMi 


-16.6 


37.4 


-1.6 


6.2 




V384 Per 


-22.3 


-20.3 


63.8 


-15.3 




AFGL 3068 


6.3 


-4.2 


11.0 


-13.3 


34.6 


TXCam 


-7.7 


56.9 


40.2 


-18.9 




GX Mon 


1.6 


-3.4 


-1.9 


3.4 




WXPsc 


-3.2 


-14.8 


23.8 


6.8 


273.7 


IKTau 


-3.3 


-1.4 


9.2 


16.4 




IRC-10529 


-3.5 


-44.8 


11.9 


38.9 





